
### Install packages

	install.packages("maptools")
	install.packages("sp")

### Load packages

	library(maptools)
	library(sp)
	library(foreign)

### Load data

	df <- read.csv("input.csv")
	attach(df)
	data <- df[,2:3]

	oh <- readShapePoly("tl_2019_36_tabblock10.shp")

### Convert to coordinates

	coordinates(data) <- ~ longitude + latitude

### Map to blocks

	res <- over(data,oh)
	write.csv(res,file="output.csv")
